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Abstract 

We construct spherically symmetric solutions to the Einstein-Euler 
equations, which contains a positive cosmological constant, say, the Einstein- 
Euler-de Sitter equations. We assume a realistic barotropic equation of 
state. Equilibria of the spherically symmetric Einstein-Euler-de Sitter 
equations are given by the Tolman-Oppenheimer-Volkoff-de Sitter equa¬ 
tion. We can construct solutions near time periodic linearized solutions 
around the equilibrium. The Cauchy problem around the equilibrium can 
be solved. This work can be considered as a trial of the generalization of 
the previous work on the problem without cosmological constants. 
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1 Introduction 

We consider the Einstein-de Sitter equation 

1 SttG 

R^V ~ ^Rah) ~ ^Ofiv = ^4 Tfin 

for the energy-momentum tensor of a perfect fluid 

= (cV -b P)U^U'' - 

Here R^u is the Ricci tensor associated with the metric 

da^ = Ofi^dx^dx'^, 

G is the gravitational constant, c the speed of light, p the mass density, P the 

pressure and is the 4-dimensional velocity. A is the cosmological constant 

which is supposed to be positive in this article. See [n§iii]- 
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Spherically symmetric solutions for the problem with A = 0 was investigated 
in [4] and the aim of this article is to describe the similar results on the problem 
with positive cosmological constants. 


In this article we suppose that the pressure P is a given function of the 
density p and pose the following 


Assumption. P is an analytic function of p > 0 such that 0 < P, 0 < 
dP/dp < c^ for p > 0, and P ^ 0 as p ^ +0. Moreover there are positive 
constants A, 7 and an analytic function fl on a neighborhood of [ 0 ,+oo[ such 
that f2(0) = 1 and 

P = Ap'>n{Ap'>-^/c^). 

We assume that 1 < 7 < 2 and - is an integer. 

7 - 1 

We are keeping in mind the equation of state for neutron stars. See a, 0 

p. 188]. 


We consider spherically symmetric metrics of the form 

ds^ = - e^^^^’^'>dr'^ - R{t, rf{d9‘^ + sin^ Odtj?). 

We suppose that the system of coordinates is co-moving, that is, 

[/O = e-^, =U^ = U^ = 0 

for = ct, x^ = r, = 9, x^ = (p. Then the equations turn out to be 

.dR 


.,-F' 


dt 
-pdp 
' dt 
dV 


= V 






-m = 


/ m 

Vp3 


+ 


47rP\ 


-R 


.,-F 


dm 

Ft 


- 1 


47r 


2Gm A 2 
— ttR 


C 2 P 


P' 


R'{p + P/c2) 


= — ttR^PV 


Here X' stands for dXjdr. 

The coefficients of the metric are given by 

P' PF'{(?pPP) = 0 


(la) 

(lb) 


(l c) 

(l d) 


and 





2Gm 

~Fr 
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In order to specify the function F, we introduce the state variable u by 

and fix the idea by putting 

with a positive constant k+ specified in the next Section. We note that there are 
analytic functions Up, Up on a neighborhood of [ 0 , +oo[ such that ^^(O) = 
flp(O) = np(0) = 1 and 

7-1 

p = Aiu^ Vt p(u ! (?) ^ 

P = AA'lu~<-^VLp{u/?). 


Here Ai := 
We put 



m = Att j pE?R'dr, 

Jo 


supposing that p is continuous at r = 0. The coordinate r can be changed to 
m, supposing that p > 0 , and the equations are reduced to 


e 


-F 


e 


-F 



m 


m 


= 1 




( 2 a) 


47rP 


c2 J 




^ + —R 
3 

p ^ -1 


dR 


?R 


c-p 


Here {d/dt)m means the differentiation with respect to t keeping m constant. 
We will change the coordinate m to r later through a fixed equilibrium, and we 
shall construct solutions near the equilibrium. 


2 Equilibrium 

Let us consider solutions independent of t, that is, F = F[r),F[ = H{r),p = 
p{r),V = 0, R = r. The equations are reduced to the Tolman-Oppenheimer- 
Volkoff-de Sitter equation 

dm 
dr 

dr 

V c'-^r 3 


= 47rr^p, 


= -[p + P!?)- 


47rr^ 




P\ - 


?K 


,/ 2Gm 


(3a) 

(3b) 
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This equation was analyzed in [^. Let us summarize the results. 

For arbitrary positive central density pc there exists a unique solution germ 
{m{r), P{r)), 0 < r <C 1, such that 


o p 2i 

m = + [r \ 2 r, 

P = -Pc- (Pc + ^’c/c^)(^^(Pc + 3Pc/c^) - + p]2- 

Here denotes a convergent power series of the form J2k>Q ■ 
We denote 


2Gm A 


K{r.m) := 1-5-—r" 

c^r o 


/ 47rr^ \ 

Q(r,m,P) := - 



(4a) 

(4b) 


we concentrate ourselves to solutions satisfyting K{r,m(r)) > 0. Moreover a 
solution (m(r), P(r)), 0 < r < r+, of (I5al)(l3bl) is said to be monotone-short if 
r+ < 00, dP/dr < 0 for 0 < r < r+, that is, Q{r,m{r), P{r)) > 0, and P —?> 0 
as r —>■ r_|. — 0 and if 


K+ := lim K{r,m{r)) = 1 — 
—0 


2G'm| 


A 

“3 


r 


2 

+ 


and 

c^A 

<5+ := lim Q(r,m{r), P(r)) = Gm+ -rl 

r^r+-0 3 

are positive. Here 


j-r+ 

m+ := lim m{r) =4 tt p(r)r^dr. 

r-)-r+-0 

We suppose that there is a monotone-short solution {fh{r), P{r)), 0 < r < r+, 
satisfying and fix it hereafter. 


As for sufficient conditions for the existence of monotone-short prolongations, 
see [S]. Anyway, the associated function u = u(r) turns out to be analytic on a 
neighborhood of [0, r-i.] and 

'“(’’) = (^-s - ^) + [^- 1 - - r]2 

r+K+ 

as r —)• r+ — 0. See [3 Theorem 4]. 


3 Equations for the small perturbation from the 
equilibrium 

Using the fixed equilibrium m = Tfi(r), we take the variable r given by its inverse 
function. We are going to a solution near equilibrium of the form 

R = r(l y), V = rv. 
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Here y, v are small unknowns. The equations turn out to be 



® dt 


{l + yf P 5 , , , 

-5- -v-^{rv) + 

p or 


-G(l+ !,)(- 


m 


+ 


47r \ c^A 

^P)+^(l + !/) 


.r3(l + ?/)3 c2 

( r'^v^ 2Gm ^ \2\ 

c2r(l + y) ^ +2/) ) X 

P \-i(l + y)2 dP 




pr dr 


(5a) 


(5b) 


Here m = m{r) is a given function and p, P are considered as given functions 
of r and the unknowns y, z(:= rdy/dr) as follows: 

p = p(r)(l + y)"^(l + y + z)~^, 

P = P(r)(l - r(u(r))(3y + z) - $(u(r),y, z)). 


Here 


p dP 
P~d^ 


and ^{u,y,z) is an analytic function of the form X)fco>o ki+k 2 > 2 '^^°y^^ ■ We 

shall denote such a function by [m; y, 2:]o;2 hereafter. Moreover we shall use 


P 

1 H —— — (1 + 

Gp 


+ +(r-i)(32/ + ^)+ 


c^p 


c^pj 


[w(r);y,z]o;2 


4 Analysis of the linearized equation 

Let us linearize (I5a1) (IFbl) : 



e = E 2 y''+ Eiy'+ Eoy, (6b) 

at 
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where y" = d'^y/dr'^,y' = dy/dr and 
E 2 = e-2^(p + p/c2)-ipr, 


El 


d 


~ ^^~P/d^P) + log(Prr^ 

dTrG 7-7— 

Eo = _3(r-l)P+ 


+ (-l- 3re-2ff + 3(r - l)e-2^^(l + P/c^p)-A [p + P/c'^)-^-'^+ 
V J r dr 


3e 

-A 


-1 


-2^(p + p/c2)-i-^pr+ 
r dr 
du^ 




Here X, XXX denote the evaluations along the fixed equilibrium. Putting 
Cy := -e^^{l + P/c^p){E 2 y" + Eiy' + E^y), 


we get the linearized wave equation 

d'^y 

dP 


+ £y = 0. 


We can rewrite C in the formally self-adjoint form 

r ^ddy 

£y = -T-r°'-r + Qy^ 

0 dr dr 


where 


a = e"+^- 


PTr^ 


l + P/c^p 


b = e' 


3H-F 


4 - 

r p 


1 -I- P/c^p 


Q = -e^^l + P/c'^pEo. 


It is easy to see that Q is bounded on 0 < r < r+. Therefore [U Proposition 7] 
is still valid: 


Proposition 1 The operator ‘Io,'D{‘Io) = G“(0,r+), To2/ = £y in the Hilbert 
space L‘^{{0,r+)\b{r)dr) admits the Friedrichs extension T, a self adjoint oper¬ 
ator, whose spectrum consits of simple eigenvalues Ai < A2 < • • • < Aj/ < • • • —>■ 
-1-00. 


As in [4], we introduce the new variable x instead of r defined by 


tan2 9 
1 -I- tan2 9 


with 
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Here 



Then there are positive constants Cq,Ci such that 

r = Co^/xil + [a;]i) as x —>■ 0, 
r+— r = (71(1 — x)(l + [1 — x]i) as x —>■ 1 


Using this variable, we can write the operator C as 





N \dy 


dy 


where Ti(x), Lq(x) are analytic functions on a neighborhood of [0,1], and 


N := 


27 

7-1’ 


which is supposed to be an even integer. Since the discussion is quite parallel 
to that of m §5], we omit the details. We may assume that ^+/7r = 1 without 
loss of generality, by changing the scale of t. 


5 Rewriting (I5al } ( I5bj) using jC 


Let US go back to the equations (f5a|h[5b|l . We shall use the analysis of dPfdr 
given in m (6.2)]: 


\ dP ^ dP 1 (9 

+ (1 + d.^/T)^-^PT{3y + z)+ 
rp or rp dr rp or 

P 1 dP 

rp rp dr 


where [(50], [(51] are given by [H (6.3a)], [H (6.3b)]. 
We put 


the right-hand side of ((5bl) = [R2] + [iil] -|- [lU]^, 

where 

(1 + P dv 

[R2] := --- -v{v + w) with w = r—, 

p or 

/ 1 P \ 

[W] := c^(l + y) - r^(l + y)\\ + P/cV)"' ( " ^^) ■ 

We put 

[Rl] = [i?3] -h [P4] -f [R5] + [P6] -f [R7] 
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as in [?]. (The symbols $^, 7 ^ of [1] are replaced by $,r in this article.) But 
the analysis of [W] is new: We put 

[W] = + [Wl] + [W2] + [W3] + [W4], 

[Wl] := c^y - r^il + y)^{l + P/c^p)~^( - 
= [WIL\ + [WIQ], 

[WIL\ :=cy- ir-y + (F - l)(3j, + z)-. 

IW2] := -r-yi + !/)‘(l + P/c7)-‘(l + a,i/r)~PT(3y + z). 

rp or 

[W3] := -r2(l + y)\l + p/cV)-'^[QO], 

rp 

1 dP 

[ 1 ^ 4 ] := -r^{l + y)\l + P/c^p)-^ 

rp dr 

Then it follows from (I3bp that 

[i?l] + [W]^ = [i?3L] + [mQ] + [R4L] + [R4Q] + [i?5] + [i? 6 ] + [R7] + 

(j 

+ {[WIL] + [WIQ] + [W2] + [W3] + [H"4])^. 

O 

Let us define 


1 + Gi = (i + a,$/r)(i + 


2Gm 




c^r(l + y) 




2Gm A gN-if + P/c^p 
c^r 3 / 1 + P/c^p 


(1 + 2 /)= 


Then we have 

-e-2^(l + P/c2p)-i£2/ = [P3L] + [P4L] + [W1L]^+ 

o 


Putting 

G 2 ■= Gi([P3L] + [R4L] + [W1L]^) + 

- ([P3Q] + [P4Q] + [P6] + [R7] + [R2])+ 
-([W1Q] + [W3] + [W4])|, 

H2 ■■= e^G2, 

Hi := 


(l + P/cV)-Hl + Gi), 












we can write 


X (the right-hand side of (fhbll i = —HiCy — H 2 . 
Using these analysis, we claim as [H Proposition 11]: 

Lemma 1 we have 


{dzHi)Cy + dzH^ = 0 mod (1 — x) 


as X ^ 1. 


Here “Qi = Qo mod (1 —x)” means that there exists an analytic function 
u}{x, y, z, V, w(:= rv'),y', y") such that Qi = Qo + (1 ~ x)u}. 

Proof is similar to that of [U Proposition 11]. We see 

{dzHi)Cy + dzH2 = e^[S] 


and we have to show [S'] = 0, where 

[S] := {dzGi)(e-^^{l + P/^p)-^Cy + [i?3L] + [i?4L] + [lUlL]A/3) + 

-h Gidz{[RiL] + \mL] + [WIL\KIZ)+ 

- dz ([i?3Q] -h [i?4Q] -h [i?6] [R7] + [R2] + {[WIQ] + [1U3] -f [lU4])A/3). 

But we have 

^ -^^(['^5] + [W2]|) - dz{[R7] + [1U4]|), 


since dz[R3L], dz[R4.L], dz[R3Q], 9,[i?4Q], a^[i?6], 9^[i?2], ^^[lUlL], dz[WlQ], 
5-,[lU3] are all = 0 clearly. By a tedious calculation, we get 


dzGi 


[[R5] + [W2]^)=dz([R7] + [W4] 


1 + Gi 

= -a2$(i + 


J. 2 y 2 


2Gm 


- r‘^{i + (1 + + z) 

3/ ro dr 


c^r(l + y) 


rp dr 


so that [S] = 0. This completes the proof. 


Putting 

J := e^(l -I- P/c^p), 

we rewrite the system of equations (l5a])(l5b)) as 

dy 


dt 

dv 


— Jv = 0, 


P HiCy P H 2 — 0 . 
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(7a) 

(7b) 







6 Main results 


Let us fix a time periodic solution of the linearized equation: 

Yi = sin{VXt + Qo)'ijj{x), 

where A is a positive eigenvalue of the operator T and ip is an associated eigen¬ 
function. We seek a solution of the form 


y = e(Yi +y), 


where 


Vi = e-^ (1 + F/c^p) 


V = £(Vi + v), 
.idYi 


dt ■ 


Then we have the equation 

= ec, with w = 
The Frechet derivative of the nonlinear operator 
D‘^{w)h = 

is given by 


( 8 ) 


FFiI 


h 

F2_ 

, with h = 

k 


Fi = ^-Jk- (^{dyj)v + {d,j)vr-^'^h 
rth 

F2 = ^+ HiCh+ 
at 


-f {{dyHi)Cy + dyH2 + {[d^Hi)Cy + d^H2)r^^h+ 

+ {{dyHi)Cy + dyH2 + dn,H2r^^k. 

Thanks to Lemma 1 and considerations as [4] we can claim that there are 
analytic functions ooi, ago, an, aio, a 2 i, 020 of x,y,dxy,d^y,v,dxV such that 


Fi = ^ - Jk + (^aoia:(l ago^h, 


d 


F 2 = 


dt 

dk 


— + FliCh + ^aiia;(l — ^ + ^a2ia^(l ~ 


dx 


Thus we can apply the Nash-Moser(-Hamilton) theorem to get 

Theorem 1 Given T > 0, there is a positive number eg such that, for |e| < eg, 
there is a solution w G (^“([OjT] x [0,1]) of (0) such that 


sup 

j-\-k<n 


f d\3 / d 

Wt) \^) ^ 


and hence a solution {y, v) of f5al 


< C{n)\£\, 

L“([0.T]x[0,l]) 

of the form y = eYi + 0{e^). 
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Note that 


R{t,r+) = r+(l + esin(V^t + 0o) + 0(£^)), 


provided that if) has been normalized as = 1 ) = 1 , and that the density 
distribution enjoys the ‘physical vacuum boundary’ condition: 


P{t,r) 


C{t){r+ — r) T-i (1 + 0{r+ — r)) (0 < r < r+) 

0 (r+ < r) 


with a smooth function C{t) of t such that 


C{t) = ( 


7-1 


Q+ 

rlK+ 


1 

7-1 


+ 0(e). 


Also we can consider the Cauchy problem 
dy 


dt 


— Jv = 0, 


dv 


HiCy + i/2 = 0, 


i=0 


= '0o(a^), V =^i(x). 


i =0 


Then we have 


Theorem 2 Given T > 0, there exits a small positive 6 such that if'ibn.ipi G 
C'°°([0>1]) satisfy 


max 

k<M 




then there exists a unique solutuon {y,v) of the Cauchy problem in (^“([OjT] x 
[0,1]). Here A is suffieiently large number. 


7 Metric in the exterior domain 

Let us consider the moving solutions constructed in the preceding section, which 
are defined on 0 < t < T, 0 < r < r_|_. We discuss on the extension of the metric 
onto the exterior vacuum region r > r_|_.Keeping in mind the Birkhoff’s theorem, 
we try to patch the Schwarzschild-de Sitter metric 

ds"^ = - — {R'^)^{d6^ + sin^ 9d<p‘^) 

from the exterior region. Here t^ = t^{t,r),R^ = R^{t,r) are smooth functions 
of 0 < t < T, r+ < r < r_|_ + (5, 5 being a small positive number, and 
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The patched metric is 


ds'^ = gooc'^dt^ + 2gQicdtdr + gndr'^ + g22{dO^ + sin^ Odcf^), 


where 


goo 


goi 


gii 


522 


= K+e 


(r < r+) 

K^dtt^f -(r+ < r) 


C K 


0 (?■ < f+) 

CK^dtt^)(drt^) ^(dtR^)(drR^) (r+ < r) 


_e 2 ff = _ 


CK^ 




2Gm A 


< 

—B? (r < r+) 

— {R^Y (r+ < r). 


-R^) {drRf 


(r < r+) 
(r+ < r) 


We require that R = R^ and drR = drR^ along r = r-|_. It is necessary 
for that (722 is of class C^. Moreover, by the same way as O Supplementary 
Remark 4], we see that 


dt^ dH'^ d'^R^ 

dt' dr ^ dr^ ’ 

at r = r+ + 0 are uniquely determined so that are of class across r = r+. 
By a tedious calculation we have 


52 


dj.2 


r++0 


d^R _ .(dRy 

dr^ r+-o~'^\dr) ’ 


where 

A = - — ((9IY± _ ^7? I ^ 1 (. ,Yl - 

(? VVc^i?^ 3 (?■ dt ) \ c^R 

Since 

V 3 / V &■ c^R 3 / r=r+-0 



r+—0 


Q+ 

d^rlnl 


Y 0 , 


we see that d'^R^/dr"^ = d^R/dr^ if and only if R = 0 at r = r+, which is the 
case if the solution under consideration is an equilibrium. 


Appendix 
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Let us describe the abstract theorem we have used. This has been established 
by the author through m, la, a, therefore the proof is not repeated here. 

First of all we introduce the following classes of functions: Let us denote by 
2 l([0,1]) the set of all functions defined and analytic on a neighborhood of the 
interval [0,1], by 2lq([0,1], [0]^) the set of all functions / defined and analytic 
on a neighborhood of [ 0 , 1 ] x { 0 } x • • • x { 0 } G such that 

,yp) = ak,-kp{x)yi^ 

fclH- \-kp>q 


The set of equations under consideration is 


with 


^ - Jix,y,z)v = 0, 
dv 

— + Hi{x,y,z,v)Cy + H 2 {x,y,z,v,w) = 0 



+ Li{x)x(l - x)— + Lo{x). 
ax 


(9a) 

(9b) 


( 10 ) 


TT 1 

Here we denote z = —, w = —. We assume: 
ox ox 


(BO): The parameters No,Ni are supposed to be greater than 4. 


The coefficients J, Hi, H 2 , Li, Lq are supposed to be of class 
2to([0,l],[0]2),ao([0,l],[0]3), 2 l 2 ([ 0 ,l],[ 0 ] 4 ), 2 l([ 0 ,l]), 2 t([ 0 ,l]), respectively, and 
their domains are supposed to include x U x U,Uq x U xU x U,Uo x U x 
U X U X U, Uo, Uo, respectively, where Uq is a neighborhood of [0,1] and t/ is a 
neighborhood of 0 . 


We suppose the following assumptions: 

(Bl): 7Li(x, 0,0, 0) = J(x, 0,0) ^ and there is a constant C > 1 such that 

i < J(x,0,0) <C 

for Vx G Uq. 

(B2): We have 

d,J = 0, {d,Hi)Cy + d,H 2 = 0, d^H 2 = 0 
as X 0 and as x ^ 1. 
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Here the meaning of ‘= 0 as a; —>■ xq’ is as follows: 

Let us denote by 2l([a;o] x [0]^) the set of all functions defined and analytic 
on a neighborhood of (a;o,0, ••• ,0) G A function / in 2l([a;o] x [0]^’) is 

said to satisfy / = 0 as a; —>■ a;o iff /(a;o, yi-, - ■ ■ , Vp) = 0 for V?/i, • • • ,yp, that is, 
there is a function fl in 2l([a;o] x [0]^’) such that 

f{x, yi,--- , yp) = {x- xo)fl(x, 2 / 1 , •• • , yp). 

In the assumption (B2), the functions under consideration are regarded as 
functions of x, y, Dy, D^y, v, Dv. Here and hereafter D stands for djdx. 


Let us fix T > 0 arbitrarily, and fix functions y*,v* G C°“([0, T] x [0,1]) such 
that all y*,z* = dy*/dx, v*, w* = dv* fdx are confined in [/ for 0 < Vt < T. We 
seek a solution y,v G C'°°([ 0 ,r] x [ 0 , 1 ]) of (I5a])(l9bl) of the form 


y = y* + y, v = v*+v, 


( 11 ) 


which satisfies 

y|t=o = 0, '0|t=o = 0. (12) 

The conclusion is: There is a small positive number S{T) and a large number 
^ such that, if 

nu« \\d(d^{y*,v*)\\L-- < 6(T), (13) 

j+k<SK 

then there exists a solution {y,v) of i9a\)0E\)ill\)U^) . 


In fact the equations for w = {y, v)'^ turns out to be 


where 


— - Jv - {AJ)v* = Cl, (14a) 

dv 

— + HiCy + {AHi)Cy* + AH^ = C 2 , (14b) 


J 

AJ 

Cl 

Hi 

AHi 

AH2 


C2 


uV 

J{x,y* + y,z* + z) with ^ 

J{x,y* +y,z* + z)- J{x,y*,z*), 
dy* 


dt 


+ J{x,y* ,z*)v*, 


Hi{x,y* +y,z* +z,v* +v,w* + w) 
dv 


with w = 


dx’ 


Hi{x, y* +y,z* + z, v* + v) - Hi{x, y*,z*,v*), 
H 2 {x, y* +y,z* + z, v* + i), w* + w) + 
-H 2 {x,y*,z*,v*,w*), 
dv^ 

—^ - Hi{x,y*,z*,v*)Cy* - H 2 {x,y*, z*,v*,w* 


(15a) 

(15b) 

(15c) 

(15d) 

(15e) 

(151) 

(15g) 
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We write the equations (I14ap(ll4bl) as *P(w) = c, where c = (ci,C 2 )^. The 
domain of the nonlinear mapping fp is it, the set of all functions w = {y, v)’^ G (Sg 
such that 


\y\ + \Dy\ + |u| + \Dv\ < eg. 


(16) 


Here eg is small so that (HU) implies y,z,v,w G U, and £ = C: x 2:, C: = 
(7°“([0,T] X [0,1]), Co = ^ I ^lt=o = V'|t=o = 0}- The Frechet deriva¬ 

tive Htp of ip is 


D^{w)h 


'Fi' 


h 


for h = 

-^2 

k 


(17) 


where 


f)h / f) \ 

Fi = — - Jk - (^{dyJ)v + idzJ)v—jh, (18a) 

f)k 

F2 = ^ + H,Ch + 
at 

-\-[(dyFli)Cy + dyFl2 + ({dzHi)Cy + dzF[2)-^'^ h + 

+ (^{dyHi)Cy + dyH2 + dwH2 ■ (18b) 


Thanks to (B2) we see that there are a^i, G 2t([0,1], [0]®),/r = 0,1,2, v = 1, 0, 
such that 

dh 

Fi = — - Jk + (agix(l - x)D + agg)h, (19a) 

F2 = ^ + H,Cy + 
at 

-|-(aiia;(l — x)D -|- aio)h + (a 2 ix(l — x)D + a 2 g)k, (19b) 


where are analytic functions of x,y, Dy, D^y,v, Dv. Under this situation, 
we can apply the Nash-Moser(-Hamilton) theorem. 
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